Statistical Mechanics of Splay Flexoelectricity in Nematic Liquid Crystals 
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We develop a lattice model for the splay flexoelectric effect in nematic liquid crystals. In this 
model, each lattice site has a spin representing the local molecular orientation, and the interac- 
tion between neighboring spins represents pear-shaped molecules with shape polarity. We perform 
Monte Carlo simulations and mean-field calculations to find the behavior as a function of inter- 
action parameters, temperature, and applied electric field. The resulting phase diagram has three 
phases: isotropic, nematic, and polar. In the nematic phase, there is a large splay flexoelectric effect, 
which diverges as the system approaches the transition to the polar phase. These results show that 
flexoelectricity is a statistical phenomenon associated with the onset of polar order. 



I. INTRODUCTION 

Flexoelectricity is a coupling between elastic deforma- 
tion and electrostatic polarization in a liquid crystalline 
medium. In general, a splay or bend deformation of the 
nematic director leads to an electrostatic polarization, 
which can be observed as a macroscopic dipole moment 
of the liquid crystal. Conversely, an applied electric field 
induces an electrostatic polarization, which leads to a 
combination of splay and bend distortions in the nematic 
director. Since its discovery in 1969 by Meyer the 
flexoelectric effect has drawn great interest because of 
its possible applications @,|| in strain gauges, transduc- 
ers, actuators, micro power generator and electro-optical 
devices. 

There have been many experimental and theoretical 
studies to determine the flexoelectric coefficients of ne- 
matic liquid crystals 0, 1, 1 0, & I, M El El El, M 
El, El, E3, EH, using a range of different approaches. 
For typical calamitic (rod-shaped) liquid crystals, the 
splay and bend flexoelectric coefficients are in the range 
of 3-20_pC/m. However, in recent experiments, Harden 
et al. [J Q found that bent-core liquid crystals have a 
surprisingly large bend flexoelectric coefficient, up to 35 
nC/m, roughly three orders of magnitude larger than the 
typical value. With this large bend flexoelectric coeffi- 
cient, bent-core liquid crystals may be practical materials 
for converting mechanical into electrical energy. 

For theoretical physics, a key question is how to ex- 
plain the large flexoelectric effect found in bent-core ne- 
matic liquid crystals, so that it can be exploited for tech- 
nological applications. Our conjecture is that the large 
flexoelectric effect is a statistical phenomenon associated 
with nearby polar phase. Near a polar phase, a nematic 
liquid crystal is on the verge of developing spontaneous 
polar order, and hence any deformation of the director 
should induce a large polar response. To test this con- 
jecture, we would like to build a model with nematic and 
polar phases, and determine the behavior of flexoelectric 
effect as a function of temperature above the nematic- 
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polar transition. In this paper, we begin the study by 
investigating the splay flexoelectric effect in a system of 
uniaxial pear-shaped molecules. In a subsequent paper, 
we will investigate the more complex case of bend flex- 
oelectricity in bent-core liquid crystals, as in the experi- 
ments. 

To study the splay flexoelectric effect, we generalize 
the Lebwohl-Lasher lattice model of nematic liquid crys- 
tals In the original Lebwohl-Lasher model, each 
lattice site i has a spin n^, which represents the local 
nematic director, with the symmetry n.; — > — n.j. In 
our generalization, the spins represent the orientations of 
pear-shaped molecules, which do not have that symme- 
try. For that reason, the interaction between neighboring 
spins includes three terms — one term favoring nematic 
order, another term favoring polar order, and a third 
term that couples polar order with splay of the nematic 
director. With this interaction, we find a phase diagram 
with isotropic, nematic, and polar phases, as illustrated 
in the snapshots of Fig. [TJ The nematic phase has a flex- 
oelectric effect, which increases as the system approaches 
the polar phase. Thus, this calculation demonstrates ex- 
plicitly that the flexoelectric effect is a collective, statis- 
tical phenomenon, which is strongest near the transition 
to a phase with spontaneous polar order. 

The plan of this paper is as follows. In Sec. II, we set 
up the theoretical framework leading to the model inter- 
action and discuss the relevant order parameters. In Sec. 
Ill, we present the Monte Carlo simulation methods and 
results, for both the phase diagram and the flexoelec- 
tric effect. In Sec. IV we present a mean-field theory for 
this model, and compare the simulation results with the 
mean-field approximation. Finally, in Sec. V we discuss 
and summarize the conclusions of this study. 



II. MODEL 

In this study, our goal is to simulate the splay flexoelec- 
tric effect in a system of uniaxial pear-shaped molecules. 
For these simulations, we construct a lattice model that 
can represent both nematic and polar order. In this 
model, the local molecular orientation at lattice site i 
is represented by a unit vector hi . If the system has ne- 
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FIG. 1: (Color online) Snapshots of the simulation results in the three phases: (a) Isotropic, (b) Nematic. (c) Polar. The 
software V_Sim is used [2p| ]. and the color of each molecule represents the polar angle 8 away from the z-axis. 



matic order, then the molecular orientations tend to be 
aligned along a preferred axis; i.e. there is a nonzero or- 
der parameter (/^(^i 'd)), where d is the overall director 
and P2 is the second Legendre polynomial. If the system 
has polar order, then the molecular orientations tend to 
point in a particular direction; i.e. there is a nonzero or- 
der parameter (Pi(n.j -d)), where P\ is the first Legendre 
polynomial. Note that the system can have nematic or- 
der without polar order, but it cannot have polar order 
without nematic order. 

The lattice Hamiltonian must have four terms: one 
term that favors nematic order, one term that favors po- 
lar order, one term that gives a coupling between polar 
order and an applied electric field, and a final term that 
gives a coupling between polar order and splay of the ne- 
matic director. The term favoring nematic order can be 
written simply as —A(hi ■ hj) 2 , summed over all pairs 
of neighboring sites i and j, as in the Lebwohl-Lasher 
model [19]. The term favoring polar order can be writ- 
ten even more simply as —B(hi -hj), again summed over 
all pairs of neighboring sites i and j, as in the Heisenberg 
model of magnetism. The coupling between polar order 
and an applied electric field can be written as — E ■ hi, 
summed over i. 

The coupling between polar order and nematic splay is 
somewhat more subtle. For this coupling we need a lat- 
tice expression for the local splay between neighboring 
sites i and j. Our expression for the local splay should 
depend only on the nematic director, and hence it should 
be invariant under the transformation h — > — ft. We can- 
not describe splay by the scalar V • h, because it is not 
invariant under that transformation. Rather, we must 
describe splay by the vector n(V • h), which has the cor- 
rect symmetry. In the following calculation, we let Latin 
letters refer to lattice sites and Greek letters refer to di- 
rections. 

On a continuum basis, the splay vector n( V • n) can be 
written in terms of the local nematic order tensor Q a fs(v), 



or equivalently in terms of the dyad n a (r)np(r), as 
-(n (3 n 7 )9 ) 3(n ct n 7 )] . 



(1) 



Hence, a lattice approximation to the splay vector be- 
tween sites i and j can be written as 



(2) 
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where f* 



Ti is the unit vector from site i to j on 



the lattice. After some algebra, this expression simplifies 
to 

[n(V • h)]ij = i [fij{rij ■ hj) - h^fij ■ hi) 

+h l (n l ■ hj)(rij ■ hj) (3) 
-hj(hi ■ hj)(fij ■ hi)]. 

Note that this expression is invariant under the transfor- 
mations hi — > —hi, hj — > — hj, and i <-> j. 

Now that we have found an expression for the local 
splay vector, we can couple it with the local polar or- 
der. The coupling term in the lattice Hamiltonian can 
be written as the dot product of the splay between sites 
i and j with the average polar order on these sites, 



V int = -C[n(V • n)]* 



(4) 



Simplifying with the use of Eq. ([3]), the coupling term 
between splay and polar order in the Hamiltonian is 

' 1 + h ■ h \ 2 
Mnt = -C ( nj ) ■ (hj - hi) . (5) 
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Combining all these terms, our final expression for the 
lattice Hamiltonian is 



H 



E 



A(hi ■ hj) 2 + B{hi ■ hj 



(6) 



-C I 1 + n » n i | . {n . ,, 



At this point, we want to use this lattice Hamiltonian 
to calculate the nematic order parameter (P2), the po- 
lar order parameter (Pi), and the average splay vector 
(n(V • n)) as functions of the parameters A, B, and C 
and the electric field E. In the following sections, we will 
do this calculation through Monte Carlo simulations and 
mean-field theory. 



III. MONTE CARLO SIMULATION 

As a first step in exploring this model, we carry 
out Monte Carlo simulations of a system of pear-like 
molecules interacting with the lattice Hamiltonian of 
Eq. ©. In these simulations, we use a simple cubic lat- 
tice of size 20 x 20 x 20. When an electric field is applied, 
it is in the z direction, so that the molecules tend to 
align along z, with splay in the x and y directions. The 
lattice has periodic boundary conditions in z, but free 
boundaries in x and y, so that it can form splay in those 
directions. 

The usual Metropolis algorithm was used for lattice 
updates. In each Monte Carlo step, a lattice site is cho- 
sen randomly, its orientation is changed slightly, and the 
change in energy AE is calculated. If AE < the move 
is accepted, and if AE > the move is accepted with 
probability exp(— AE/ksT). Starting from the high- 
temperature isotropic phase, the system is cooled down 
slowly with temperature steps of AT = 0.02. The final 
configuration at each temperature is taken as the initial 
configuration for the next lower temperature. Typical 
runs take about 10 5 steps to come to equilibrium, while 
runs near phase transitions take about 6 x 10 5 steps. The 
nematic and polar order parameters and the splay vector 
are calculated and time-averaged during the production 
cycle. 

The nematic order parameter (P2) is calculated by the 
usual method using the 3D nematic order tensor 
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where a and (3 = x,y,z, and N is the total number of 
lattice sites. The largest eigenvalue of this order tensor 
corresponds to (P2). 

To calculate the polar order parameter (Pi), we as- 
sume that polar order is oriented along the same axis as 
nematic order, as expected for uniaxial molecules. The 
eigenvector corresponding to the largest eigenvalue of the 



nematic order tensor Q a p is the instantaneous director 
d. Hence, the polar order parameter is calculated as the 
average dot product of the director with the molecular 
orientation, 



(Pi) 



1 N 
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(8) 
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The splay vector is calculated from Eq. (|5J), averaged 
over the four bonds in the (x, y) plane. The magnitude of 
this vector gives the average angle between the molecular 
orientations on neighboring lattice sites. For that reason, 
we report this magnitude as (A9). 

Figure [5] shows plots of the order parameters (P2}, 
(Pi), and (Ad) as functions of temperature for several 
values of the interaction parameters. In Fig. [Ha), for a 
small polar coupling B and no applied electric field, we 
see an isotropic-nematic transition at high temperature 
followed by a nematic-polar transition at low tempera- 
ture. At the isotropic-nematic transition, the nematic 
order parameter goes from zero to a nonzero value. Here 
the transition is rounded by finite-size effects; we would 
expect a sharp first-order transition for an infinite sys- 
tem. Throughout the nematic temperature range, the 
polar order parameter and splay are both zero. At the 
nematic-polar transition, the polar order parameter be- 
comes nonzero, and this polar order induces an accom- 
panying splay. The nematic order parameter decreases 
as the system moves into the polar phase, because the 
splayed molecular orientation partially averages out the 
alignment, as shown in the snapshot of Fig. [IJc). 

In Fig. [2fb), for a larger polar coupling B, we see a 
direct transition from the isotropic to the polar phase, 
with no intervening nematic phase. In this case, the ne- 
matic and polar order parameters both become nonzero 
at the same transition temperature. Once again, the po- 
lar order induces a splay, which inhibits the growth of 
the nematic order parameter. 

The simulation results for the phase diagram are shown 
in Fig. [DJc). In this phase diagram, the vertical axis 
shows temperature while the horizontal axis shows the 
polar coupling B for a constant nematic coupling A. For 
small B the phase diagram shows isotropic, nematic, and 
polar phases, with a nematic range that decreases as B 
increases. At a sufficiently large value of B, the nematic 
phase disappears and there is a direct transition from 
isotropic to polar. Note that this phase diagram is quite 
similar to the phase diagram found in recent work on the 
2D isotropic, tetratic, and nematic phases [2lj . 

The question is now: What happens to the nematic 
phase when an electric field is applied? To answer this 
question, Fig. Hfd) shows the simulation results for the 
same parameters as Fig. (Ha), but in the presence of a 
small electric field. In the high-temperature isotropic 
phase, the field induces some polar and nematic order, 
but this effect is very small. However, in the nematic 
phase, the field induces a more substantial polar order, 
and that polar order induces a splay in the nematic di- 
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FIG. 2: Monte Carlo simulation results for the order parameters (Pi), {P2), and (A#) as functions of temperature T, for 
different values of the interaction parameters chosen to show various types of transitions: (a) Zero-field results for A — 1.5, 
B = 0.09, C — 0.3, showing the isotropic-nematic and nematic-polar transitions, (b) Zero-field results for A = 1.5, B — 0.4, 
C — 0.3, showing the direct isotropic-polar transition, (c) Phase diagram for zero field, (d) Simulation with applied electric 
field E — 0.06, for the same parameters as in part (a), showing the induced polar order and splay in the nematic phase. 



rector, i.e. a converse flexoelectric effect. Both the po- 
lar order and the splay are quite temperature-dependent, 
increasing as the system approaches the nematic-polar 
transition temperature, as would be expected for a diver- 
gent susceptibility above a second-order transition. The 
nematic-polar transition is rounded off by the applied 
field, and the polar order parameter and splay saturate 
in the low-temperature polar phase. 

To provide further insight into the effect of an applied 
electric field, Fig.[3]shows the splay as a function of tem- 
perature for several values of the field. Within the ne- 
matic phase, the splay increases as the electric field in- 
creases, as expected for the converse flexoelectric effect. 
This trend is reasonable because an increasing electric 
field enhances polar order. For small field the splay is 
quite sensitive to temperature, but for large field it be- 
comes less temperature-dependent, as the induced polar 
order grows larger and approaches saturation. In the 
low-temperature polar phase, the splay shows the op- 
posite trend with electric field; it now decreases as the 
field increases. This trend is reasonable because an in- 



creasing electric field cannot enhance the polar order, 
which is already saturated; it only aligns the direction of 
polar order. This alignment reduces the induced splay, 
since splay necessarily involves some misalignment of the 
molecular orientation. 

For comparison, Fig. [4] presents a plot of splay as a 
function of temperature for several values of the inter- 
action coefficient C. This graph shows that the splay 
increases as C increases, over the full temperature range, 
in all phases. This result is reasonable because the co- 
efficient C represents the flexoelectric coupling between 
polar order and splay. 

As a final point, note that the behavior presented here 
can only occur in the limit of small splay A6> < ir/N, 
where N is the system size. In the opposite limit 
A8 > tt/N, the system is too large for one single splay 
from side to side. Instead, it must break up into modu- 
lated structures consisting of regions of splay separated 
by domain walls. These modulated structures might 
be splay stripes or even more complex two- or three- 
dimensional arrangements of splay cells *22j. We have 
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FIG. 3: Variation of splay as a function of temperature for 
several values of the applied electric field. 




FIG. 4: Variation of splay as a function of temperature for 
several values of the interaction coefficient C, with a small 
field E = 0.06. 



observed such modulated structures in our simulations, 
but we have not explored them in detail because they are 
not likely to occur in experiments, where the magnitude 
of splay is generally small. 



A. Ising model 

For a simple Ising-type model of the splay flexoelec- 
tric effect, we suppose that the system has well-defined 
nematic order, with variable amounts of splay and po- 
lar order. Consider a particular site i surrounded by 
six nearest neighbors on a cubic lattice. We suppose 
that site i has its director along the z-axis, as do the 
two neighbors above and below, while the four near- 
est neighbors in the xy-plane have directors that are 
splayed outward by a small angle A9. The polar or- 
der at any site j is represented by an Ising spin vari- 
able <7j = ±1, which indicates whether the molecular 
orientation is pointing up or down along the local di- 
rector. Thus, the central site i has the molecular ori- 
entation fii = (7^(0,0, 1), while the six neighbors have 
the orientations n+ x = er+ x (sin AO, 0, cos AO), n_ x = 
ct_ x (— sin A9, 0, cos AO), n +y = cr +y (0,sinA#, cosA#), 
n_ y = <7_ y (0, — sin A#,cos A#), n +z = er+ z (0, 0,1), and 
n-z =<t_ z (0,0,1). 

We now substitute these expressions for the molecular 
orientations into the lattice Hamiltonian of Eq. ([5]). As 
usual in mean-field theory, we assume that all the neigh- 
bors of site i have polar order given by (crj) = M. (This 
quantity is called Pi in the other sections; here we use 
the symbol M to emphasize the analogy with the Ising 
magnetization.) The mean potential experienced by site 
i, expanded to second order in the small splay A9, is then 



-A(6-4(A<9) 2 ) - BMai(6 
-2CA0{a l + M) - Eu t . 



2(A9) 2 



(9) 



Hence, the effective field acting on the Ising spin <jj is 

E cS = E + 6BM + 2CA9 - 2BM(A9) 2 . (10) 

As a result, the polar order parameter must satisfy the 
self-consistency equation 



M 



((jj) = tanh 



k B T 



(11) 



tanh 



E + 6BM + 2CA9 - 2BM{A9f 



IV. MEAN-FIELD CALCULATION 



In this section we discuss two approximate analytic 
approaches to solve the problem. First, we map the in- 
teraction onto an Ising model, and use this Ising model to 
calculate the splay and polar order as functions of tem- 
perature and electric field. Second, we present a more 
general mean-field calculation with full rotational degrees 
of freedom, and use it to calculate the full phase diagram 
with isotropic, nematic, and polar phases. 



Furthermore, minimization of the mean potential over 
the splay A9 gives 



A9 = 



CM 



2A + BM 2 ' 



(12) 



Solving Eqs. (fTTj) and (fT2"|) simultaneously gives the equi- 
librium values of the splay A9 and polar order M, as 
functions of electric field E, temperature T, and ener- 
getic parameters A, B, and C. 

To calculate the response to an electric field in the 
nematic phase, we assume that E, M, and AO are all 
small, and expand Eqs. (| 1 1 [) and (|12p to linear order in 
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this quantities. These equations imply 



M = 
AO = 



E 



k B T- (6B + C 2 /A) 7 

CE 

2A[k B T-(6B + C 2 /A)Y 



(13) 
(14) 



Note that Eq. (T3"|) gives the polar order parameter in- 
duced by an applied electric field, while Eq. (|14[) gives 
the converse flexoelectric effect induced by the field. Both 
of these responses increase as the temperature decreases 
toward the second-order nematic-polar transition at the 
temperature 



C 2 

kBT NP = 6B + — 



(15) 



At this transition, they diverge as (T — T/vp)~ 7 , with 
critical exponent 7 = 1, as expected for the susceptibil- 
ity to an applied field, in mean-field theory for the Ising 
model. 

For a more precise calculation, we solve Eqs. (fTTj) 
and (11211 numerically as functions of temperature and 
field. The numerical results for splay A9 and polar order 
M are shown in Fig. [5] As in the approximate analytic 
calculation above, we see a second-order nematic-polar 
transition. The high-temperature nematic phase has no 
polar order or splay without a field, but an applied field 
induces both of these quantities. By contrast, the low- 
temperature polar phase has both spontaneous polar or- 
der and spontaneous splay, and they both increase mod- 
erately when a field is applied. 

Although the Ising model is successful in explaining 
some features of our Monte Carlo simulations, it is in- 
complete because it assumes perfect nematic order — the 
molecules can have only two possible orientations, up and 
down. It cannot describe the behavior of the nematic or- 
der parameter as a function of temperature. For that 
reason, we proceed to a more general mean-field theory, 
in which each molecule has full rotational degrees of free- 
dom. 



B. General Mean-Field Calculation 



In mean-field theory, the free energy can be written as 



F = U-TS={H} + k B T {log p), 



(16) 



averaged over the single-particle distribution function p. 
Thus, our goal is to express the single-particle distribu- 
tion function in terms of some variational parameters, 
calculate the energetic and entropic terms in the free en- 
ergy, and then minimize the free energy over those vari- 
ational parameters. 

For this mean-field calculation, we write the molecular 
orientation at each lattice site in terms of the polar angle 
Oi and azimuthal angle <pi with respect to the local direc- 
tor. We assume the distribution function depends only 



on the polar angle 9i, and hence write 

exp[«iPi(cos6'i) + v 2 P 2 (cos Oi 



JJ'exp^iP^cos^) + V2P2 (cos 0i)]d0i 



(17) 



where v\ and v 2 are variational parameters. The or- 
der parameters are then (Pi) = J P 1 (cosO)p(0)dO and 
(P 2 ) = Jq P2(cos 0)p(9)d9, and the partition function is 
Z = J™ exp[t>iPi(cos#) + V2P2 (cos 9)]d9. The entropic 
contribution to the free energy is therefore 

-TS = k B T{\og P {9 t )) = k B T[v l (P l )+v 2 (P 2 )-\og(Z)] 

(18) 

per lattice site. 

As in the previous section, we suppose that site i has 
its director along the z-axis, as do the two neighbors 
above and below, while the four nearest neighbors in the 
xy-plane have directors that are splayed outward by a 
small angle A9. To calculate the average energy, we com- 
bine our distribution function of Eq. (|17p with the lattice 
Hamiltonian of Eq. ©. After averaging over all the an- 
gles, neglecting terms involving the third-order Legendre 
polynomials, and expanding to second order in the small 
splay AO, we obtain 



(H) 



-A-A(P 2 ) 2 [2-2(A6>) 2 
-P(Pi) 2 [3 - (A0) 2 ] 

-c ( pW 5 + 1( ™ 



(19) 
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- P(Pi) 



1 - 



(AOf 
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per lattice site. 

We now have an expression for the total free energy 

F = -A - A(P 2 ) 2 [2 - 2(A0) 2 \ - B(Pi) 2 [3 - (AO) 2 ] 



-C(Pi)A9 



16(P 2 ) 



15 



- P(Pi) 



1 - 



(AO) 2 



+k B T[v 1 (P 1 ) +v 2 (P 2 ) -log(Z)} 



(20) 



per lattice site. In this expression, note that (Pi), (P2), 
and Z are all determined by the parameters v\ and V2 
in the distribution function. Hence, the free energy is a 
function of just three variational parameters: v±, V2, and 
AO. Thus, in the mean-field calculation, we must min- 
imize the free energy numerically with respect to those 
three parameters. After this minimization, we can calcu- 
late the order parameters (Pi) and (P 2 ), and hence de- 
termine whether the system is in an isotropic, nematic, 
or polar phase. 

Figure [6] shows the numerical mean-field results for the 
order parameters (Pi) and (P 2 ) and splay AO, as well as a 
complete phase diagram as a function of temperature T. 
For a small polar interaction B — 0.09 there are two tran- 
sitions, first from the high-temperature isotropic phase 
((Pi) = 0, (P 2 ) = 0) to the intermediate nematic phase 
(((Pi) = 0, (P 2 ) 7^ 0), and then from the nematic phase 
to the low-temperature polar phase ((Pi) ^ 0, (P 2 ) ^ 0). 
The isotropic-nematic transition is first-order, while the 
nematic-polar transition is second-order. On increasing 
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FIG. 6: Numerical mean-field calculations of the order parameters (Pi), (P2), and A8 as functions of temperature T, for 
different values of the interaction parameters chosen to show various types of transitions: (a) Zero-field results for A — 1.5, 
B = 0.09, C — 0.3, showing the isotropic-nematic and nematic-polar transitions, (b) Zero-field results for A = 1.5, B — 0.4, 
C — 0.3, showing the direct isotropic-polar transition, (c) Phase diagram for zero field, (d) Simulation with applied electric 
field E — 0.06, for the same parameters as in part (a), showing the induced polar order and splay in the nematic phase. 
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the polar interaction strength B, the polar phase becomes 
stable even at higher temperature. For B = 0.36, these 
two transitions merge into a single first-order transition 
directly from the high-temperature isotropic phase to the 
low-temperature polar phase. If E ^ 0, the polarization 
and splay are nonzero even in the nematic phase and 
scale with the magnitude of the field. For nonzero field, 
the magnitude of the splay increases on reducing temper- 
ature and is enhanced greatly near the transition to the 
polar phase. 

Note that the numerical mean-field results of Fig. [5] 
are very similar to the Monte Carlo simulation results 
of Fig. [2l both in the overall phase diagram and in the 
splay and polar response to an electric field. This simi- 
larity demonstrates that the mean-field theory captures 
the essential physics of this model. 

V. DISCUSSION 

In conclusion, we have developed a lattice model for 
splay flexoelectricity in a system of uniaxial pear-shaped 
molecules. This model predicts a phase diagram showing 
isotropic, nematic, and polar phases, and it further pre- 
dicts a converse flexoelectric effect in the nematic phase. 
The converse flexoelectric effect is proportional to the 
applied electric field, and it increases dramatically as the 
temperature decreases toward the nematic-polar transi- 
tion. Indeed, we can regard this effect as a susceptibility 
to an applied field, which diverges at the second-order 
transition to a polar phase. Thus, flexoelectricity is not 
just a molecular effect arising from the microscopic inter- 



action of liquid crystals with a field. Rather, it is a sta- 
tistical effect associated with the response of correlated 
volumes of molecules, which increases as one approaches 
a polar phase. 

The recent experiments of Harden et al. @, _f] have 
found an anomalously large bend (rather than splay) flex- 
oelectric effect in systems of bent-core liquid crystals. 
We speculate that the same considerations discussed in 
this paper can explain the large bend flexoelectric coef- 
ficient in those experiments. The bent-core liquid crys- 
tal should be close to a polar phase, with order in the 
transverse dipole moments of the molecules. As a result, 
there should be large correlated volumes of molecules, 
leading to a high susceptibility to an applied field, which 
induces both polar order and bend. The detailed theo- 
retical model for bend flexoelectricity will necessarily be 
more complex than the model for splay flexoelectricity 
presented here, because the bent-core molecules are not 
uniaxial and hence their orientations must be character- 
ized by two vectors or three Euler angles. This model 
for bend flexoelectricity will be the subject of a future 
paper. 
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